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Abstract: 



We study the impact of the "intrinsic" hadron transverse momentum on the pre- 
asymptotic behavior of the diffractive electroproduction of longitudinally polarized 
p-meson. Surprisingly, we find the onset of the asymptotic regime in this problem to 
be rather low, Q 2 ~ 10 GeV 2 where power corrections due to the transverse momen- 
tum do not exceed 20 % in the amplitude. This drastically contrasts with exclusive 
amplitudes where the asymptotics starts at much higher Q 2 = 50 — 100 GeV 2 . The 
sources of such unexpected behavior are traced back to some general (the quark- 
hadron duality) as well as more silent (properties of higher dimensional vacuum 
condensates) features of QCD. 
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1 Introduction 



Over the last couple of years, following the paper H, the diffractive electroproduc- 
tion of the light vector mesons has been investigated within the framework of per- 
turbative QCD. Since this paper a few more important results have been obtained. 
First of all it is a rigorous QCD-based proof @ that an amplitude of diffractive elec- 
troproduction of vector mesons has required factorization properties. The second 
important step in the same direction is an introduction of a new type of functions 
(the so-called antisymmetric structure functions) || which are relevant objects for 
such kind of problems. Therefore, the standard perturbative QCD (pQCD) approach 
to this new class of problems which are currently intensively studied experimentally 
0] has obtained a solid basis for the future investigation. 

We should note that the diffractive electroproduction had been discussed many 
times before the recent paper [Q], see e.g. p| and references therein. However, 
all previous approaches to the problem were mainly based on some kind of quark 
models with inevitable to such methods specific assumptions not following from 
QCD. This is certainly a self-consistent approach for a heavy quark system (where 
the size of a hadron is of the order of 1/mq and is always small by kinematical 
reasons), but not for the light hadrons, where the problem of the separation of small 
and large distance physics cannot even be formulated in an appropriate way within 
the quark model. Let us remind that the main idea of the QCD-based approach is 
just such a separation of the large and small distance physics, the so-called Wilson 
Operator Product Expansion (OPE). At small distances one can use the standard 
perturbative expansion due to the asymptotic freedom and smallness of the coupling 
constant. All nontrivial, large distance physics is hidden in nonperturbative hadron 
wave functions (WF's) in this approach. They cannot be found by the perturbative 
technique, but rather should be extracted from elsewhere. Therefore, the problem 
is reduced to the analysis of the bound states within OPE. 

The problem of bound states in the relativistic quantum field theory with large 
coupling constant is, in general, an extremely difficult problem. Understanding the 
structure of the bound state is a very ambitious goal which assumes the solution 
of a whole spectrum of tightly connected problems, such as confinement, chiral 
symmetry breaking phenomenon, and many others which are greatly important 
in the low energy region. Fortunately, a great deal of information can be obtained 
even in absence of such a detailed knowledge. This happens in high energy processes 
where the only needed nonperturbative input are the so-called light-cone WF's with 
a minimal number of constituents. In this case a problem becomes tractable within 
existing nonperturbative methods. 

We have therefore a well-formulated problem of the diffractive electroproduc- 
tion of vector mesons. The formulation is based on the solid background of QCD. 
However, calculations of the pQCD approach refer, strictly speaking, only to asymp- 



1 



totically high energies |TJ]. Therefore, as usual, the main question remains: at 
what energies do the asymptotic formulae of Ref. QXJ] start to work? The prime goal 
of the present work is to attempt to answer this question. We should note that 
there are two very different types of pre-asymptotic corrections to the asymptoti- 
cally leading formula First, there exist the so-called shadowing (rescattering) 
corrections || which we do not discuss in this paper at all. Our concern is rather 
with a second kind of corrections which are related to the internal hadron structure 
of the vector meson. This kind of corrections was discussed for the first time in 
Ref. 0, and we will discuss the relation of the paper @ to our results later on. 

However, before deepening in details, we would like to recall that a very similar 
question about an applicability of the QCD-based approach to exclusive amplitudes 
has been discussed during the last fifteen years . The history of this development is 
very instructive and we believe it is worthwhile to mention it shortly here. 



As is known, at asymptotically high energies the parametrically leading contribu- 
tions to hard exclusive processes can be expressed in terms of the so-called distribu- 
tion amplitude f8|, which itself can be expressed as an integral / d 2 k±ip(k 2 L , x) 



with nonperturbative wave function ip(kj_,x), see reviews [fJ,|T(J f° r details. Distri- 
bution amplitudes 4>(x) depend only on longitudinal variables Xi and not on trans- 
verse k\ ones. The same is true for inclusive reactions where structure functions 
depend on x, but not on fc^Q. Thus, any dependence on k\ gives some power correc- 
tions to the leading terms. Naively one may expect that these corrections should be 
small enough already in the few GeV 2 region. These expectations are mainly based 
not on a theoretical analysis but rather on the phenomenological observation that 
the dimensional counting rules, proposed in early seventies (see [TTJ ), agree well 
with the experimental data such as the pion and nucleon form factors, large angle 
elastic scattering cross sections and so on. This agreement can be interpreted as a 
strong argument that power corrections are small in the few GeV 2 region. 



However, in mid eighties the applicability of the approach || at experimentally 
accessible momentum transfers was questioned |T^|,[|T^]. In these papers it was 
demonstrated, that the perturbative, asymptotically leading contribution is much 
smaller than the nonleading "soft" one. Similar conclusion, supporting this result, 





lq] , where the direct 
3GeV 2 . This method 



came from the different side, from the QCD sum rules 
calculation of the form factor has been presented at Q 2 < 
has been extended later for larger Q 2 < lOGeV 2 fl6||,[1I7| with the same qualitative 
result: the soft contribution is more important in this intermediate region than the 
leading one. 



* The formal reason for that can be seen from the following arguments. At large energies 
the quark and antiquark are produced at small distances z ~ l/Q — ► 0, where Q is typical 
large momentum transfer. Thus one can neglect the z 2 dependence everywhere and one should 
concentrate on the one variable zQ ~ 1 which is of the order of one. One can convince oneself that 
the standard Bjorken variable x is nothing but Fourier conjugated to zQ. 
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Therefore, nowadays it is commonly accepted that the asymptotically leading 
contribution to the exclusive amplitudes cannot provide the experimentally observ- 
able absolute values at accessible momentum transfers. If we go along with this 
proposition, then the natural question arises: How can one explain the very good 
agreement between the experimental data and dimensional counting rules if the 
asymptotically leading contribution cannot explain the data for experimentally ac- 
cessible energies? A possible answer was suggested recently |18| and can be formu- 
lated in the following way: very unusual properties of the transverse momentum 
distribution of a hadron lead to the mimicry of the dimensional counting rules 
by the soft mechanism at the extended range of intermediate momentum transfers. 
Numerically, the soft term is still more important than the asymptotically leading 
contribution at rather high Q 2 ~ 50 4- IWGeV 2 . 



Now we come back to our original problem of the hard diffractive electroproduc- 
tion. Having the experience with exclusive processes in mind, one could expect a 
similar behavior (i.e. a very slow approach to the asymptotically leading prediction) 
for the diffractive electroproduction as well. The main goal of the present work is 
to argue that this naive expectation is wrong. The asymptotically leading formula 
starts to work already in the region Q 2 ~ 10 GeV 2 , where power corrections do not 
exceed the 20% level in the amplitude. This surprising result is a consequence of 
very special properties of the WF in the k\ variable (as we mentioned earlier, the 
k\ dependence determines corrections to the leading term) . At the same time, as we 
argue below, these properties of WF are very sensitive to the QCD vacuum struc- 
ture, more exactly to some higher dimensional vacuum condensates. As a result, 
we find that the pre-asymptotic behavior of the p-meson diffractive electroproduc- 
tion amplitude essentially depends on the (non-) factorizability of particular mixed 
quark-gluon condensates. In our opinion, this makes diffractive processes an excel- 
lent laboratory for testing our current understanding of the QCD vacuum structure. 
In contrast to the standard QCD sum rules for static properties of hadrons Plfl 
which probe most crude characteristics of the QCD vacuum like (ipip) > (9 2 G 2 ), 
diffractive amplitudes are remarkably influenced by far more silent vacuum features. 
Moreover, in the latter case all these theoretical ideas can in principle be given a 
dynamical test by varying the photon virtuality Q 2 and measuring a deviation from 
the asymptotic predictions of Ref.fl]] (provided, of course, all other sources of the 
Q 2 dependence such as shadowing corrections |3j and the Q 2 evolution of the gluon 
distribution are properly taken into account). One can therefore say that by study- 
ing the pre-asymptotic behavior of diffractive amplitudes we actually study the fine 
structure of the QCD vacuum. 



Our presentation is organized as follows. In Sect. 2 we collect relevant definitions 
and discuss general constraints on the nonperturbative p-meson WF stemming from 
general principles of QCD. Sect. 3 deals with the calculation of the lowest moments 
(k 2 ^) and (k\) of the p- meson transverse momentum distribution from the equations 
of motion and QCD sum rules. The main observation of this section is a presence of 
strong fluctuations in the transverse direction (k\)/(k 2 L ) 2 ^> 1. We show how this 
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ratio is fixed by particular mixed quark-gluon condensates. In Sect. 4 we propose a 
model WF satisfying the general constraints of Sect. 2 as well as the lowest moments 
(k\) and {kjj of Sect. 3. All these results are applied to the study of pre-asymptotic 
effects due to the "intrinsic" transverse momentum (the Fermi motion) in the hard 
diffractive electroproduction in Sect. 5, where we also present our understanding of a 
fundamental difference between exclusive and diffractive processes. As a by-product 
of our study, we propose a simple reformulation of the diffractive amplitude in 
terms of a current Green function. Appendix A establishes a connection between the 
transverse momentum distribution we are dealing with and higher twist distribution 
amplitudes of the standard ligh-cone OPE approach. Methods for estimates of 
relevant higher dimensional vacuum condensates are described in Appendix B. Sect. 6 
presents a summary of our results. 



2 General constraints on the nonperturbative wave 
function 

The aim of this section is to provide necessary definitions and establish some essen- 
tially model independent constraints on the non-perturbative WF \I/(£, kj_) which 
follow from the use of such general methods as dispersion relations, duality and 
large order perturbation theory. For definiteness, we will consider the longitudinally 
polarized charged p-meson. 

Let us start with the definition of the non-perturbative WF \I/(£,& 2 ) in the so- 
called impact parameter representation as an infinite series of local gauge invariant 
operators (a low normalization point /i 2 ~ 1 GeV 2 is implied and b 2 = — z 2 ) 

]T^(0|J(0) 7 ^ V„) n «(0)|p(<7)> = f P m p e^(zq, b 2 ) 

= f P m p e^ J\e~^(£,b 2 ) (1) 

Here V v =Vu — Vu , Vv=du —igA v and stands for the polarization vector. 

In the asymptotic limit of hard exclusive processes [f| |9| the quark and antiquark 
are produced at small distances z ~ 1/Q — > where Q is a typical large momentum 
transfer. In this case the only remaining dependence on the variable zq ~ 1 is given 
by the leading twist distribution amplitude (DA) 0(£) : 

(0|d(0) 7 ^ V„)M0)|p(g)> = f p m p el x \zq) n (C} = f p m p e^{zq) n £^£"0(0 

f = 2u - 1 , z 2 = (2) 

In the infinite momentum frame (IMF) q z — > oo the DA 0(£) describes the dis- 
tribution of the total longitudinal momentum q z between the quark and antiquark 
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carrying the momenta uq z and (1 — u)q z , respectively. In what follows we will use 
the both variables £ and u , u = 1 — u interchangeously. 

A portion of model independent information on the DA 0(£) can be obtained 
using the dispersion relations and quark-hadron duality To this end, we study the 
asymptotic behavior of the correlation function 

i f dxe^(0\TJl(x), Jo(0)|0) = (zq) n+2 I n (q 2 ), 4 = d lv z v {i V„ z^u (3) 



(here the label || stands for the longitudinal direction which is selected by the projec- 
tor z^ : z 2 = , z^q^ ^ .) At large q 2 — > — oo the exact answer for the correlation 
function ) is given by the perturbative one-loop diagram 

I roo ds Imir\s) = 3 

7r Jo s — q l An(n + l)(n + 3) 

On the other hand, in the IMF the correlation function (|]) can be parametrized in 
terms of intermediate hadron states, the lowest one being the p— meson. Therefore 
the following relation takes place 



1 r S P 1 r°° 

- / dslml(s) p n ert = - / dslml(s) p n , (5) 

TC JO TV JO 



The only assumption made here is that the p-meson gives a non-zero contribution 
in (H), i.e. the duality interval Sp(\\) ( the subscript || specifies the longitudinal 
direction ) does not vanish for arbitrary n. No further assumptions on highly excited 
states (like those imposed in the QCD sum rules approach) are needed. In this case 
Eq. (D yields 

3^ 



-<oo / 



ffKn-oo)^^ (6) 

( we have used the fact that for the longitudinally polarized p-meson in the IMF 
= qp/m p + 0(l/q z ) .) This formula unambiguously fixes the end-point behavior 



of the DA UTD| 



(D = £ d^ n m ~ yn 2 , m - ±1) - (i - e 2 ). (7) 

We want to emphasize that the constraint (|7|) is of very general origin and follows 
directly from QCD. No numerical approximations were involved in the above deriva- 
tion. Pre-asymptotic as q 2 — > — oo perturbative and non-perturbative corrections 
are only able to change the duality interval in Eq. (^) (which is an irrelevant issue, 
anyhow) but not the parametric 1/n 2 behavior. 

The same analysis can be done for the transverse distribution. It moments are 
defined analogously to Eq.(||D through gauge invariant matrix elements 

(0\d^(i V„ t u ) 2n u\p(q)} = f pmp£ ^{-t 2 T { ^^ (k 2 r) , (8) 
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where transverse vector = (0, t, 0) is perpendicular to the hadron momentum 
% = {%■> 0j_, q z ). The factor j^jij is introduced to @ to take into account the inte- 
gration over <p angle in the transverse plane: / d(f>(cos <fi) 2n / J d<p = (2n — l)!!/(2n)!!. 
By analogy with a non-gauge theory we call (k\) in this equation the mean value of 
the quark perpendicular momentum, though it does not have a two-particle interpre- 
tation. The usefulness of this object and its relation to the higher Fock components 
will be discussed at the end of this section. Here we note that Eq. (||) is the only pos- 
sible way to define (k±) in a manner consistent with gauge invariance and operator 
product expansion. 

To find the behavior (k 2 ™) at large n we can repeat the previous duality argu- 
ments with the following result []: 



This behavior (for the case of pion) has been obtained in Ref.||19|| by the study of 
large order perturbative series for a proper correlation function. Dispersion relations 
and duality arguments transform this information into Eq.(^). It is important to 
stress that any nonperturbative wave function should respect Eq.@ in spite of the 
fact that apparently we calculate only the perturbative part ( see the comment after 
Eq.([?D). The duality turn this perturbative information into exact properties of the 
nonperturbative WF. 



The most essential feature of Eq. (|9]) is its finiteness for arbitrary n. This means 
that higher moments 

(pn) = J dk 2 ± d£k?V (10) 

do exist for any n. In this formula we introduced the nonperturbative k±) 
normalized to one |. Its moments are determined by the local matrix elements (§). 
The relation to the longitudinal distribution amplitude </>(£) looks as follows : 

J dkmz, k\) = 0(0 , £ d&(o = i (ii) 

The existence of the arbitrary high moments (k 2 ™) means that the nonperturbative 
k±) ^ ans a ^ large transverse momentum k\ faster than any power function. 
The relation fixes the asymptotic behavior of k'jj at large k\. Thus, we 
arrive at the following constraint: 

(fcf ) = J dk 2 ± d^k 2 ^{k 2 ± ,0 ~ n\ , n -> oo (12) 

' Here and in what follows we ignore any mild (nonfactorial) n-dependence. 

*To avoid possible misunderstanding, we stress that our definitions are very different from those 
made in the light-cone perturbation theory || where (fc^)ss = j\d^J d 2 k±k 2 _\ li i>Bs(£,,k'j_)\ 2 
(here BS stands for bound state). Our approach does not suppose at all the existence of a bound 
state equation in QCD, but rather relies entirely on the logic of OPE. 
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We can now repeat our duality arguments again for an arbitrary number of 
transverse derivatives and large (n — > oo) number of longitudinal derivatives. The 



result reads EC 



dk 2 jf^(ki, e - ±1) ~ (i - e) k+i (is) 



The constraint (|13"D is extremely important and implies that the k\ dependence 
of ty(kj_,£) comes exclusively in the combination k\/{l — £ 2 ) at £ — > ±1. This 
means that the standard assumption on factorizability ^/(k"j_,^) = i)(k±)<f>(£) is at 
variance with very general properties of the theory such as duality and dispersion 
relations. The only form of k±) satisfying all the constraints (|7|) , (|T2|) and 
is the Gaussian with a very particular argument : 



^{k\ -+ oo, ~ exp - A? . x — (14) 




(here A 2 is a mass scale which can be fixed by calculating the moments (kj_), (k\) 
etc.) Strictly speaking, so far we have only established the validity of Eq.(|i~4"D in a 
vicinity of the end-point region £ —» ±1. In Appendix A we discuss in what sense 
Eq.(|i~4"|) can be approximately valid in the whole range of the £ variable. 



Here a few comments are in order. Our result (J14Q essentially supports the old 
quark model-inspired idea on the SU(6) invariance of hadron WF's. Indeed, the very 
same analysis carries without any changes for pseudoscalar 7r - and K -mesons (non- 
zero s-quark mass effects are subleading). A possible violation of SU(6) is expected 
via a variation of the parameter A 2 for different hadrons. Moreover, the ansatz (|14"D 
coincides with the well known harmonic oscillator WF, except for the mass term. 
One should stress that there is no room for such term in QCD. Its inclusion violates 
the duality constraint (|7|) since in this case we would have 

^ exp(— y/n) , n — ► oo (15) 

instead of the 1/n 2 behavior (|7J). In other words, a true nonperturbative WF must 
respect the asymptotic freedom which is incompatible with the quark model -type 
mass term in the WF |. This difference leads to interesting phenomenological con- 




sequences, see Ref.[]T8| for more detail. Furthermore, we would like to clarify the 
correspondence between our results ( |14"1) and those obtained in Ref. |23| on similar 
grounds. There, it was also claimed that k\ and £ enter the WF only in the combi- 
nation fcl/(l-£ 2 ) but a different WF was suggested, k]_) ~ Q{s Q -k 2 J '(l-£ 2 )). 



§To set this more accurately, one can say that a possible (scale-dependent) mass term in k 



must renormalize to zero at a normalization point /i 2 ~ a few GeV 2 where the duality arguments 
apply. This conclusion is not at variance with popular models for the QCD vacuum such as e.g. 
the instanton vacuum (see p2[ for review). 
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Our approach differs from that of Ref. [2~3| in that we apply the duality arguments 
for currents with derivatives and in addition use the large order perturbation theory 
analysis. As a result, we end up with the Gaussian instead of the G -function. 

Our final comment concerns with the meaning of Eq.(^). As has been explained 
above, (k\) is not literally an interquark transverse momentum. The notion of 
interquark transverse distance is senseless in a gauge theory because quark transverse 
degrees of freedom are undistinguishable from longitudinal degrees of freedom of 



higher Fock states at any normalization point due to equations of motion, see Ref. |24 
and Appendix A. Thus an attempt to separate them would be at variance with 
exact equations of motion. Any theoretically consistent calculation of higher twist 
effects must include these two effects simultaneously at each given power of Q 2 . 
However, the two contributions may well differ numerically. In fact, in all known 
examples where higher twist effects for exclusive amplitudes were calculated |]| [T| 
a contribution coming from the ipGip Fock state turns out smaller than a two-particle 
contribution due to the transverse motion by typically a factor of 2-3 , and comes 
with an opposite sign. It is probably worth noting that in the light-cone QCD sum 
rule approach to the pion form factor at intermediate Q 2 a two-particle DA of twist 
4 yields the same 1/Q 4 behavior as the leading twist DA, while a contribution due 



to TpGip DA is down by the extra power of Q 2 |T7]]. Thus, one can expect that a 
calculation based on the use of the light-cone WF k±) alone can serve as a 
reasonable estimate of higher twist effects. It is just the line of reasoning taken 
in this paper. In no sense we claim that a WF like (|T4]) exhausts the p-meson 
properties, nor we pretend to give a complete evaluation of higher twist effects in 
diffractive electroproduction (see Sect. 5). 



Lowest moments (k±) , (k\) and vacuum con- 
densates 



The general constraints of the previous section are insufficient for building up a 
realistic nonperturbative WF for the p-meson. To fix the latter, we follow the same 
logic as in the analysis of the distribution amplitudes and calculate the lowests (kjj 
and (k\) moments of k±). The physical meaning of the second moment (k'jj 
is clear: this quantity serves as a common scale for power corrections in physical 
amplitudes, cf. the next section. The fourth moment (k±) produces a next-to- 
leading power correction and signals on how strongly k±) fluctuates in the k\ 
plane. In this section we will calculate the lowest moments (kj_) and (kj_) by a 
combined use of the equations of motion and QCD sum rules technique. With the 
knowledge of (kj_) and (k\) we then construct a model WF kj_) which will be 
used to estimate higher twist effects in diffractive electroproduction in Sect. 5. 
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We start with a general Lorentz decomposition of the matrix element 

\ ( MT^ V« i V/3 u\p{q)) = e { ^q a qpf p m p (u 2 ) 

perm 

+ (4W + 4"W - 2e^q a qp)A 2 (16) 

(here perm stands for the symmetrization in respect to a, (3). We have defined 
the last term such that to reproduce the usual definition of the first term as the 
second moment of the leading twist DA <p{u) when all Lorentz indices are contracted 
with the light-like vector z^. The second term gives rise to the definition (^|) after 
contraction with the transverse vector t p . With the definition (|16| ) at hand, our 
strategy is to relate the unknown parameters (k 2 ^ , Ai , A% to some matrix elements 
of gauge invariant operators. At the second step these matrix elements are estimated 
using the QCD sum rules. 



The first relation between the coefficients Ai,A 2 is obtained by multiplying 



2m p 



(ki) -5A 1 -A 2 = (17) 



Another relation can be derived when one multiplies Eq.fllB"!) by q p . This yields the 
following constraint 

(4 A) ga + 4 A) ^) m p(A + A 2 ) = -(0\d lp (gG pa V/3 +gG pP V a )u\p{q)) 

= -(ei% + 4\)m;i/ (18) 

The matrix element fllS]) parametrized by the number rj can be estimated by studying 
the following correlation function 

i J dx e iqx (0\T{u lp d(x) dj^gG^ V/3 +gG p(3 V o )«(0)}|0> (19) 

The matrix element of interest appears as the lowest intermediate hadron state 
contribution to the imaginary part of (|19|) . On the other hand, the correlation 
function fll9|) can be calculated for q 2 — > — oo in QCD. The perturbative contribution 
can be easily read off a very similar sum rule in Ref. p4| and numerically turns out 



negligible. The first power correction ~ (g 2 G 2 ) in ( |T9"D vanishes, and we end up with 
the following simple estimate 



2 



rj ~ —a s ^\ ~ 1.8 x GeV 2 (20) 
27 f p m s p 

We next contract Eq. (|i~6[) with 8 a p. Then the matrix element in the LHS becomes 
(0\d^(i Vafu\p(q)) = -(0\dgG pulul5 u\p(q)) = -e^W p \ (21) 
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The number A = 24 ± 3 MeV was calculated a long time ago [25J by the QCD sum 
rules method. Using (PD , (0), (0) we finally obtain 



(*1> 



2m p 
3/ p 



[m p A + f p m p (u 2 ) + 377] ~ (420 Mel/) 2 



(22) 



Note that the last term in Eq. (|22|) is rather small in comparison with the other two 
terms (for the second term we used here the asymptotic DA <p{u) = 6u(l — u) with 
{u 2 ) = J duu 2 (f>(u) = 0.3 ). This fact is quite understandable as the matrix element 
in Eq.flllf) is proportional to first moments of three-particle qGq DA's while their 
normalizations are of the order of m p X. Therefore, an uncertainty in Eq.(pO|) does 
not play an important role in our estimate (0). We will use this observation in the 
calculation of (k\). Numerically {k\) for the p— meson (^) is somewhat larger than 
the analogous parameter for the pion {kj_) n ~ (330 MeV) 2 ]T0| ( see also Eq.( |A.9| ) 
in Appendix A). 



The same type of analysis can be done for calculation of the fourth moment (k\) 
though the algebra in this case becomes somewhat more complicated. We write 

— ^ (0|<#V V« i V/3 % Va i Ve u\p(q)) = e^lqaqpqxqzfpmpiu 4 ) 

perm 

+(6 a p6x^ + perm)-f p rn p (kj) + {8 a pq\q^ + perm)C x ] 

o 

+ {qfi £ a ' \pq\qt + P erm - qaqpqxqi)C 2 + q^e^q^x^ + perm)C 3 (23) 



+ (S fia ef ) q x q^ + perm)C A + (5 pa 5p X e^' + perm)m 2 C 5 

Proceeding analogously to the previous case and neglecting numerically small con- 
tributions due to second moments of three-particle DA's, we can obtain a system of 
linear equations on the coefficients Cj. The solution reads 



(A) 



1M) 

24 m r 



C 2 
C 3 

c 4 
c 5 



f P m 3 p {u 4 ) 



20 

48 m p 
1 

2o /p 77 "* 
48 m 



We finally contract Eq.(p3|) with S^Sx^ and make use of 
equation holds 



(24) 



-^.{gaGfuW)) + ±{0\d lp g 2 Glpu\p(q)) 



e^m 2 p 



|). Then the following 



--fp^) + 2^(k[) 



m, 



(25) 
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We have thus reduced the calculation of the fourth moment (k\) to the calculation 
of the matrix elements of operators containing only the gluon tensor G^ v and not 
the gluon potential itself as in the original problem (^). To estimate these matrix 
elements, it proves convenient |20) to consider the following nondiagonal correlation 



functions 



T t (p) = % J dxe< px (0\ua a pd(x) d lfl g 2 G^u(0)\0) 
T 2 (p) = i [ dxe ipx (0\ua a/3 d(x) J7^w(0)|0) (26) 



Perturbative contributions vanish in the chiral limit in the both correlation functions 
(p6|) due to an odd number of gamma matrices. Thus, a leading as p 2 — ► — oo 
behavior comes from condensate terms in the OPE for Ti,T 2 . Dividing T\ by T 2 
we get rid of the p-meson residue ~ {Q\ua a pd\p(q)) in (|26|). Then for the reduced 
matrix elements Ai, A2 defined by 

(0\d^g 2 Gl p u\p(q)) = ej^mpn , (0| d^(gaG) 2 u\p(q)) = ^m 5 p X 2 (27) 

we obtain after a simple algebra the following relations for the matrix elements of 
interest in terms of mixed vacuum condensates of the dimension 7 : 



Ai 



/, (0\u(gG) 2 u\0) 

(0\uu\0) 
f p (0\u(gaG) 2 u\0) 
m\ (0|w|0) 



:28i 



Here the second of Eqs.(P^) is obtained analogously to (^). Using Eq.(P5|) we finally 

(Sd - ^ra> { u>) + i mem + im^Gim (29) 

v ±; 10 pV 1 4 <0|uu|0> 8 (0|mm|0) V ; 

This is the main result of this section: we have explicitly expressed (k\) in terms of 
the vacuum expectation values (VEV's) of the dimension 7 operators. Naively, one 
could estimate these condensates by factorizing them into the products of the quark 
(ipijj) and gluon (g 2 G 2 ) condensates. This procedure, based on the factorization 
hypothesis, does not work in the given case: there are essential deviations from 
the factorization prediction in VEV's (p8[) . The non-factoraziblity of mixed quark- 
gluon matrix elements of such type has been studied in by two independent 

methods with full agreement in estimates between them. The first one |26|] was based 
on the analysis of heavy-light quark systems which allows to obtain restrictions on 
the VEV's like (|28|), while the second method has related the vacuum condensates 
of the form (|28|) to some pion matrix elements known from PC AC. As our results 
essentially depend on the magnitude of the condensates we shortly review the 



method of Ref.pO| in Appendix B. Here we only formulate the result of this analysis. 
A measure of non-factorizability introduced by the correction factors Ki , K 2 in the 
matrix elements 

(0\u(gG) 2 u\0) = ^(01^10) <0|uu|0> 
o 
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(0\u(gaG) 2 u\Q) 



1 



K 2 (0\g 2 G 2 \0) (0\uu\0) 



(30) 



( K\ = K2 = 1 in the factorization limit) is approximately the same K\ = K% = 
K ~ 3 for the both mixed operators appearing in Eq.(|28|). A possible uncertainty 
of this estimate does not exceed 30 
following numerical estimate 



20| . Using Eqs. (p9|) and ( p0[) we arrive at the 



(*1> 



j( 9 2 G 2 ) + ^( 9 2 G 2 ) + ^m^) 



0.14 GeV 4 



(31) 



This large value of (kj_) seems rather surprising as it implies strong fluctuations of 
the nonperturbative WF \l/(£, k\). To measure the magnitude of these fluctuations 
it is convenient to introduce a dimensionless parameter 



R 



4.2 if (k\) ~ (420 MeVf 



(32) 



The phenomenon of strong fluctuations in \l/(£, k'jj has been found previously by 



analogous methods for the pion wave function in Ref. |2(J with approximately the 



same number for the ratio (32). In terms of the QCD vacuum structure these 
fluctuations are due to the numerical enhancement of the high dimensional quark 
gluon mixed condensates or, what is the same, the large magnitude of the parameter 
K in Eq.pl). 



4 Model wave function ^(w, k\) 



The results obtained so far are essentially model independent. We have fixed the 
form of the high-fc^ tail of the true nonperturbative WF \l/(£, k\) Eq. (|l4D by the use 
of the quark-hadron duality and dispersion relations. Furthermore, we calculated 
the lowest moments {k 2 } and (kj_) using the equations of motion and QCD sum 
rules. Now our purpose is to build some model for the true nonperturbative WF 
k±) which would respect all general constraints of Sect. 2 and incorporate the 
effect of strong fluctuations in the transverse k\ plane found in the previous Sect. 3 
(see Eq.(^)). In spite of the fact that this inevitably brings some model dependence 
into any physical quantity calculated within such a wave function, there remains a 
region where the results of such a calculation become practically model independent. 
This happens when the answer is essentially determined only by the lowest (k'jj and 
(k\) but not higher (kj_), etc. moments of the WF. As we shall see in Sect. 5, this 
is precisely the case in the hard diffractive electroproduction for sufficiently high 
Q 2 > 10 GeV 2 . 

The effect of strong k\ fluctuations fl32|) cannot be taken into account by a simple 
Gaussian ansatz of the type ( |14D which would yield the factor of 15/7 in Eq.(^). 
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The only possible way to satisfy Eq.(|3"2"D is to have a second pick in k 2 ) which 
can have a small magnitude but reside far away from the the first Gaussian one. To 
reconcile the existence of this second pick with the asymptotic behavior ( fH] ) it must 
enter as a pre-asymptotic term which would provide the constraint (B3). Following 



Ref.|T8| we suggest the following ansatz for \l/(£, k\ 



^0, k]_) = A I exp 



i(3 2 uu 



cexp 



^(3 2 uu 



If 



(33) 



In this ansatz the parameters c and I determine the magnitude and position of the 
second pick, respectively. We have found that the best fit corresponds to the choice 
(3 = 125 MeV , c = 0.15 , Z = 30 . For these values of parameters R ~ 4.1 , (kj) ~ 
0.13. For u ~ 1/2 the second pick is located at k\ ~ 1 GeV 2 . The existence of 
two characteristic scales in \l/(£, kj_) ( 0.4 GeV and 1 GeV) is rather interesting and 
may probably be linked to phenomenological models for the QCD vacuum like the 
instanton vacuum model []22] and/or the constituent quark model. If one does not 
insist on (or disbelieves) the effect of amplified fluctuations in the transverse k\ 
plane (|32|), then the second exponent in Eq.(p3|) should be omitted. It is important, 
however, that in order to have the same second moment {k\) the parameter (3 must 
be simultaneously rescaled to (3 ~ 330 MeV . In this case the WF takes the following 
form : 



^(w, k 2 ± ) = A exp 



k\ 



\(3 2 uu 



(3 = 330 MeV 



(34) 



We will use in Sect. 5 the both forms in order to reveal the importance of large (k\) 



(|31"1) and, respectively, of the non-factorizability of VEV's (p0|) in the pre-asymptotic 
behavior of the diffractive amplitude. 



The model WF (^) meets all general and numerical constraints derived in pre- 
vious sections of this paper. The integration of (|3~3D over k\ yields the asymptotic 
form of the distribution amplitude <p(u) = 6u(l — u) which is known to be a good 
approximation to the true non-perturbative DA of the p-meson. Furthermore, nu- 
merical coefficients we have obtained for the ansatz (BBf) are not very different from 



those found in |L8| for the case of the pion WF. Our conclusion is that the transverse 
momentum distributions in pion and p-meson are to a large extent alike. 



5 Hard diffractive electroproduction 



In this section we come back to our initial aim. The results of proceeding sec- 
tions expressed in condensed form by Eq.(^) will be used for the study of the 
pre-asymptotic effect due to the Fermi motion in diffractive electroproduction of 
the p-meson. Our prime goal is to get an estimate for the onset of the asymptotic 
regime in this problem. 
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The applicability of perturbative QCD (pQCD) to the asymptotic limit of the 
hard diffractive electroproduction of vector meson was established in Ref.[|l| using 
the light-cone perturbation theory. The authors of [|I| have proved that for the 
production of longitudinally polarized vector mesons by longitudinally polarized 
virtual photons the cross section can be consistently calculated in pQCD. There 
was found that at high Q 2 the amplitude factorizes in a product of DA's of the 
vector meson and virtual photon, the light-cone gluon distribution function of a 
target, and a perturbatively calculable on-shell scattering amplitude of a qq pair off 
the gluon field of the target. After the factorization of gluons from the target the 
problem seems to be tractable within OPE-like methods. 



It would be desirable to have an explicit realization to the idea of treating the 
process of interaction of the virtual photon with the gluons from the target and a 
sequent conversion into a vector meson as some kind of a current Green function, 
to which methods of the light-cone OPE could be applied. Such reformulation has 
proved to simplify enormously the analysis of the pion form factor and some other 
exclusive amplitudes. This approach is particularly convenient for the study of 
higher twist effects in exclusive processes, where calculations of light-cone pQCD 
become tremendously complicated. Here we propose such a method. Though no 
new result will be given, we believe it is still worthwhile to present it in this paper. 
The reason for this is several-fold. First, our method reproduces the results of Ref.|l]] 
more easily and moreover extents them by taking into account an asymmetry of the 
gluon distribution of the target. Second, it allows to establish a connection with 
the nonperturbative p-meson WF defined within the OPE. Note in this respect that 
the question on what WF should be substituted in the amplitude is of peculiar 
interest and causes some controversy in the literature (see below on this). Third, 
the approach we suggest is ideally suited for a complete evaluation of twist 4 effects 
including three-particle qGq DA's (see discussion in Sect. 2). Though we do not carry 
such a calculation, it can be done within this formalism provided the qGq DA's are 
known. Thus, our reformulation will be used below merely to rederive the formulae 
of Ref.0, within OPE. This will enable to show that the WF which enters Eq.(g]]) 
below is in fact the soft light-cone WF (10) defined as a set of matrix elements which 
has been discussed in the previous sections of this paper. 



Let us start with the general form of the amplitude as a matrix element of the 
electromagnetic current = e(2/3-U7 M w — l/Sd^^d) : 

M = e fl (N(p — r)p(q + r)\j /1 (q)\N(p)) , (35) 

where e M is the polarization vector of the photon (only the longitudinal polarization 
is considered, see [[!]]) and r stands for the momentum transfer. We will neglect 
the masses of the nucleon and p-meson in comparison with the photon virtuality 
Q 2 : m 2 p = m 2 N = 0. For the momentum transfer r we consider the limit r 2 = 0, 
but r M ^ 0. We next note the following. The factorization of the gluon field of 
the target [l], @ means that that these gluons act as an external field on highly 
virtual quarks produced by the photon with Q 2 — > oo. This statement can be 
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given a formal meaning by introducing into the QCD Lagrangian an additional 
term SS = ig J dzj't{z)G c t where j° = qy^t a q is the color current and stands for 
the external gluon field from the target. Expanding the path integral to the second 
order in this external field corresponds to the two-gluon exchange : 



M=-faJ dzdy(N(p-r)\(%(z)<%(y)\N(p)^ 

(36) 

One has to stress that our Eq.(^) is not another "proof of factorization, but rather 
should be considered as a reformulation of results of Ref. @] (and also ||, see Eq.([39|) 
below) in a way which is very suitable for practical calculations. Following Ref.0, 
we find it convenient to choose the light-cone gauge for the external gluons 

rod 

G^z) = n v / daG^(z + an) , n^G^ = (37) 

J 

in which the nucleon matrix element can be parametrized in terms of the asymmetric 
gluon distribution function || : 

I at( \\r<a( \nbi \\at( \\ ^ u(p - r)hu(p) ( p a n p +p p n a \ 



7V2-1 (-4)(pra) V M 
Fi(X) 



Here one exploits the fact that in the limit r 2 = the momentum transfer is pro- 
portional to the nucleon momentum r M = £p^ with £ = Q 2 /C2pq) being the Bjorken 
variable. Under this circumstance the fractions of the total nucleon momentum car- 
ried by the gluon are xp^ + yr^ = (x + £y)p^ = Xp^ and (x - £y)p^ = (X - 
where y = 1 — y . It is thus convenient || to parametrize the gluon distribution 
in terms of the asymmetric gluon distribution F^(X) which goes over to the usual 
gluon distribution function Xf g (X) in the symmetric limit £ — > 0. 

Using the translational invariance and the fact that u(p — r)nu(p) / (—Apn) = 
\/l — £ , we finally present the amplitude in the form 

M = -N[^l{ 9 ^ (^j^-tJo dX X(X-Z + i e) 

dzdye+* (+e^ x -^y + c -««) (p(q + r)|T{j M (z)^(0)^(y)}||0> (39) 



We have thus reduced the amplitude of diffract ive electroproduction to the three- 
point correlation functions of quark currents between the vacuum and the p-meson 
states. One can see that this correlation function is dominated by the ligh-cone and 
thus a leading contribution is given by a combination of non-local quark -antiquark 
string operators of twist 2 whose matrix elements yield the leading twist DA of the 
p-meson. Coefficients in front of these operators are due to most singular parts of 
quark propagators near the light-cone. Retaining only this contribution to (B5|), we 
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recover the asymptotic answer of Ref.JIJ in the form presented without derivation 



in 



4irV2ea sfp 1 f ^WVR f 1 <f>(u) 
M ~ N c Q Jo dX X(X-£ + i 5) Jo 1* (40) 

where <f>(u) is the standard light-cone DA (0) (redefined to the interval u = [0, 1] ) 
whose asymptotic form is <fi as {u) = 6uu. Remind here that it is defined as <p{u) = 
J dk\^(u, k\) where the soft WF ty(u, k\) is introduced by Eq.(|l~0|). 



Going over to higher twist corrections to the asymptotic answer (f40[ ), one has to 
distinguish between two classes of contributions. Corrections of the first class come 
from keeping less singular terms in the quark propagators in (|39"1). Relevant terms 
in the light-cone expansion for the quark propagator contain an additional gluon 
which is to be attributed to the final meson and expressed in terms of three-particle 
DA's. A corresponding calculation is technically involved, though it is feasible within 
the formalism presented here. The second class of corrections is due to the quark 
transverse degrees of freedom (the Fermi motion). Calculating only this contribution 
we make an educated guess on the scale of higher twist corrections in the diffractive 
electroproduction (see the comment at the end of Sect. 2). Technically, retaining this 
correction amounts to the multiplication of (|40|) by the following correction factor 
(following the notations of Ref. |7j we reserve the symbol T(Q 2 ) for the correction 
in the cross section) 



.H&I? d^{u,k]) (n2 J H ^ 2 1-2 



k 2 ^ /(uu) 

T (D2\ - n4__ - (Q'+ffV(uS))2 V~ -{Q*+P ± /(uu) y 

where ^>(u,k"j_) was defined earlier in terms of the local matrix elements (|T0|). By 
definition, in the asymptotics Q 2 — > oo VT = 1. Deviations from \/T = 1 determine 
a region of applicability of the asymptotic formula P0| ) . 

Now we are in position to discuss numerical estimates for the correction factor 
(PI) in order to answer the main question formulated in Introduction. The authors 
of have observed that the choice of ^f(u, k±) in a factorized form $f(u, A^j = 
0(ii)-?/>(/?[) leads to a very slowly raising function JT(Q 2 ) which approaches 0.8 at 



rather high Q 2 > 40 GeV 2 depending on the model chosen for ^f(u, k\). On the 
other hand, as we stated earlier, Eq. (|39"D unambiguously demonstrates that it is the 
soft nonperturbative WF (|T^) that has to be substituted in Eq. fl4"T|) . Note in this 
respect that any perturbative tails like l/k 2 ^ in the WF are by definition absent 
within the logic based on a use of OPE. Such contributions effectively reappear only 
as radiative corrections to leading 0{a° s ) results. It is clear that these corrections 
are always small unless there are special reasons which make them parametrically 
leading as Q 2 —>■ oo (this is the case e.g. for the pion form factor). Using the model 
WF ( P^D we obtain a very fast approach to the asymptotics in (f4"l~p. The correction 

factor y / T(Q 2 ) reaches the value 0.8 already at Q 2 ~ 10 GeV 2 that implies a rather 
low onset for the asymptotic regime in the amplitude (|39|). 
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Figure 1: The correction factor T(Q 2 ) in the cross section for different choices of 
the p-meson WF ^(u, k±). 




The technical reason for this behavior is quite clear: large values of k\/{uu) 
are exponentially suppressed by the WF (p3|) and thus this term remains much 
smaller that Q 2 for Q 2 > 10 GeV 2 . The result of our calculation is shown for the 
correction T(Q 2 ) in the cross section on the Figure (the thick curve). To demonstrate 
the importance of amplified fluctuations in the k\ — plane ([32]) (and thus of the 
non-factorization property (30) ), we also show (the thin curve) the result of the 



calculation when the WF ( |3~4"D not including the second pick is used. One sees 
that the difference between the two curves disappears for Q 2 > 40 GeV 2 , while 
it becomes quite sizeable (20 - 30 % ) in the region Q 2 = 5 — 10 GeV 2 . One may 
therefore speculate that the non-factorizability of higher dimensional condensates 
(|3T)|) could be tested experimentally by measuring the cross section in this region. 
Still, one has to note that prior to any quantitative prediction, a more refined analysis 
is required in this region, including a complete evaluation of twist 4 corrections due to 
the (ipGip) Fock state. For sake of comparison with the discussion of Ref. J7J we also 
plot (the dashed curve) the factor T{Q 2 ) corresponding to the choice |7] of the WF in 
the factorized "dipole" form ^(w, k 2 ± ) ~ uu/(k\ + /j 2 ) 2 for fi = 300 MeV. Note that 
for such WF the second moment (k±), as defined by Eq.flTUD, does not exist at all , 
and can only be introduced by means of the bound state equation, see the footnote 
after Eq.flTD]). The choice \i = 300 MeV corresponds to (k 2 ± ) BS ~ (220 MeV) 2 (this 
number follows from nowhere and is taken exclusively for illustrative purposes). 
What matters, however, is the fact that all other options discussed in @ for a 
factorized form of WF ^(u, k\) = 4>(u)i()(k\) ( including the above ansatz with 
{k 2 ± ) B s - (300 MeV) 2 and (k 2 ± ) BS ~ (600 MeV) 2 and some other choices) all lie 
under this curve. One can therefore see quite unambiguously the difference in the 
correction T(Q 2 ) for the factorized ^(u, kj_) = <f)(u)il)(k\) and the non-factorized WF 
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in the whole range of Q 2 . Furthermore, this difference is critical for the issue of 
the onset of asymptotic regime which is of our prime interest in this paper. While for 
WF ( |3"3"D the correction \JT in the amplitude reaches the value 0.8 already at Q 2 ~ 
10 GeV 2 (or equivalently T{W GeV 2 ) ~ 0.64), for the "dipole" ansatz #(u, ~ 
uu/{k\ + /i 2 ) 2 with /i = 300 MeV this happens only at Q 2 ~ 35 Ge^ 2 . For other 
choices of a factorized WF ^(u, k\ ) = 0(m)-?/'(/c 2 _) discussed in a corresponding 
value is pushed further to even higher Q 2 > 40 GeV 2 . On the other hand, at 
Q 2 = 10 GeV 2 the difference between our prediction for T(Q 2 ) and models discussed 
in J7| constitutes at least a factor of two. 



Our use of the ansatz (^) can, however, meet objections as apparently it yields 
a large contribution to the integral ([II]) coming from the region of small k\ which 
is expected to be cut off on physical grounds. Here we remind that the ansatz (|33"D 
has been suggested to satisfy the general constraints of Sect. 2 and the particular 
values of the moments (Sect. 3). On the other hand, exclusive amplitudes at large 
Q 2 are generally only sensitive to the high-fc^ tail of the WF which dominates the 
moments (&j_)> (k±) etc. Thus one expects that these are only global characteristics 
like (k"j_) that matter in Eq. pi]) . It is easy to see that this is indeed the case as soon 
as the WF is defined to possess all moments, see Eq.(TO). Powers of l/(uu) arising 



when the integrand in (flT]) is expanded in a series over k\ are harmless due to the 



general constraint ( |T4"D which excludes end-point singularities of the u-integration. 
In this case the u— and k\— integrations in Eq . fl4i"D decouple with the exponential 
in Q 2 (in practice this means Q 2 > 5 GeV 2 ) accuracy, and the answer is therefore 
directly expressed in terms of moments 

(£l> nm<*a> 1 3.4 GeV 2 28.4 GeV A 



Tm = 1 -^ + 2W ^ + ... = 1 -—— + —— + ... (42) 

Therefore our approach is self-consistent in the sense that any question on a large 
distance contribution to the integral fl4"T| ) is lifted up to the definitions of the mo- 
ments, as has to be expected on general grounds. Here we emphasize again that 
our definitions of moments are specific to OPE, see the footnote after Eq . (|i~0|) . The 
numbers in Eq.(^) illustrate the observed fast growth of T(Q 2 ): power corrections 
tend to cancel each other before they become small separately. On the other hand, 
it can be seen from Eq.(|42|) that for Q 2 < 10 GeV 2 higher moments of kj_) turn 
out essential and thus in this region the answer for T(Q 2 ) is to some extent model 
dependent. 



Thus our final conclusion is that the onset of the asymptotic regime for diffractive 
electroproduction of the longitudinally polarized p-meson is approximately Q 2 ~ 
10 GeV 2 where corrections due to the quark transverse degree of freedom constitute 
less than 20 % in the amplitude. This can be traced back to the fact that in the case 
at hand the power corrections are given by the matrix elements of local operators 
and in fact are fixed completely by the independent calculation of the moments. This 
is the consequence of the structure of Eq.(|41"|) and the fact that the WF is a function 
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of the single variable k\j (uu). Only global, but not local characteristics of ^(u, k\) 
in the k\ plane are important. This situation can be confronted with the case of 
exclusive processes. In the most well studied problem of the pion form factor the 
asymptotic regime has been found to be pushed further to Q 2 ^> 10 GeV 2 [17 



where the subleading "soft" contribution is still larger than the leading asymptotic 



one. One reason for this [17j is that a special power correction comes due to the 



Feynman end-point mechanism for the form factor which provides the whole answer 
~ 1/Q A in absence of radiative corrections. The asymptotic answer ct s (Q 2 )/Q 2 is 
in fact the radiative correction to this soft term which nevertheless wins in the 
asymptotics Q 2 — > oo. In our opinion, this is precisely where the most fundamental 
difference between exclusive and diffractive amplitudes lies. The academic Q 2 — > oo 
limit of the former is given by the hard rescattering mechanism which is 0(a s ) and 
thus is superseded at all available energies by the nonperturbative soft contribution 
which is down by the power 1/Q 2 , but at the same time does not contain a smallness 
due to a s . On contrary, in the diffractive amplitude power corrections which are 
again 0(a° s ) ( but not attached to the end-points of the u-integration) compete with 
the asymptotic amplitude which itself is also 0{a° s ). It follows from this argument 
that the onsets of the asymptotic regime for the exclusive and diffractive amplitudes 
differ parametrically by the factor l/a s . This conclusion is not in contradiction 
with the fact that for the pion form factor the nonleading soft contribution is 
larger than the asymptotic answer up to very high Q 2 = 50 — 100 GeV 2 . We are thus 
convinced that the analogy with the pion form factor problem does not work in the 
case at hand. There are no reasons to expect the onsets of the asymptotic regime to 
be similar in the diffractive electroproduction and exclusive processes. On contrary, 
they differ parametrically in l/a s . Moreover, the explicit calculations suggest that 
the asymptotic regime in the p-meson diffractive electroproduction starts already 
at Q 2 ~ 10 GeV 2 . We stress that this conclusion refers only to the diffractive 
electroproduction of the longitudinally polarized p-meson, the situation with the 
transverse polarization or diffractive charmonium production can be quite different. 



6 Summary 

We hope that we have succeeded in explaining a few facts by this paper. Below we 
collect them together. 

1. Very general principles of the theory such as the dispersion relations and quark- 
hadron duality allows one to fix the asymptotic form of the light-cone WF k 2 ^). 
It is the Gaussian of the single argument ^(w, k\) ~ exp(— k^/S^uu) , while the 
mass term in WF or/and the factorized WF k±) = </>(/u)^>(&jJ are absolutely 
forbidden by the above principles. 

2. The light-cone WF k±) does not have a two-particle interpretation and is 
related by equations of motion to higher twist distribution amplitudes. However, 
there exists a numerical difference between their contributions, which allows to re- 
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tain only the former to estimate an importance of higher twist effects for a particular 
process. 

3. Properties of ^(u, k±) are strongly influenced by the fine structure of the QCD 
vacuum. More precisely, the p-meson WF ty(u,k±) is very sensitive to the non- 
factorization of quark gluon vacuum condensates of dimension 7. 

4. The pre-asymptotic behavior of the p-meson diffractive electroproduction ampli- 
tude is critically depending on properties of the WF ty(u, k^_). When the correct 
WF is substituted in the amplitude, the asymptotic regime starts at rather low 
Q 2 ~ 10 GeV 2 , in drastic contrast to exclusive processes where the onset of the 
asymptotic regime is pushed to much higher Q 2 = 50 — 100 GeV 2 . 

5. This vast difference between exclusive and diffractive processes is not accidental, 
but rather has the clear origin : the asymptotic regime for the former is delayed in 
comparison to the latter parametrically by the factor l/a s (Q 2 ). 

6. In the intermediate energy region Q 2 = 5 — 10 GeV 2 the diffractive amplitude 
is sensitive to the non-factorization of dimension 7 vacuum condensates. Thus, the 
fine structure of the QCD vacuum is probed there. By preliminary estimates, the 
non-factorization has an effect of the order of 20 - 30 % in the cross section. 

7. To get a quantitative description of the diffractive amplitude for intermediate en- 
ergies Q 2 = 5 — 10 GeV 2 a more refined calculation is required, including a careful 
analysis of a model dependence of predictions for such Q 2 and complete evaluation 
of twist 4 effects with account for three-particle qGq distribution amplitudes. 

8. The latter technically involved calculation can be considerably simplified using 
our representation of the diffractive amplitude in terms of the current Green func- 
tion, Eq.(|39|) which re-expresses the results of 0, [3] in a form convenient for practical 
purposes. Necessary for this calculation qGq distribution amplitudes for p-meson 



will be discussed elsewhere 128 
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Appendix A 



The purpose of this appendix is to establish the connection of the problem of trans- 
verse momentum distribution with the standard light-cone operator product ex- 
pansion (OPE) and explain how (kj_) can be directly expressed in terms of some 
three-particle matrix elements. We consider this problem on the example of the 
pion. This choice is motivated by the fact that an analysis of higher twist distribu- 
tion amplitudes in this case is considerably simplified in comparison to the p-meson 
due to smaller number of independent Lorentz structures. The complete set of 
three-particle DA's of twist 4 beyond the asymptotic regime has been constructed 
in Ref.[^] using the (one-loop) conformal symmetry and the QCD sum rules tech- 
nique. Here we use these results to show that the combination b 2 uu ( or k\/{uu) 
after the Fourier transform ) is a single argument in the pion WF to the leading 
conformal spin accuracy. This is a justification for our choice of the model WF fl33|) 
where the only argument k\/uu appears. 



The twist 4 three-particle DA's of interest are introduced as follows 



{®\d(-x)~i a ~i b gG^(vx)u(x)\ir(q)) = q a {q^x u - q v x^) — / D(aiqx)$\\(ai) (A.l) 

qx 1 



+ 



qx 



av ' ) 

qx 



D(aiqx)§j_(ai) 



Here 



J D( ai qx) = J da 1 da 2 da 3 5(l -a l -a 2 - a 3 ) e - iqxiai - a2+va3) (A.2) 
To the leading conformal spin accuracy (i.e. in the asymptotic form) the DA's 



$\\(ai),$±(ati) are |24j 



lO/n-5 2 0(3(0(1 — 02) 





where 5 2 is defined through the matrix element 

(0\dj v igG vli u\ir(q)) = fj 2 q^ , 5 2 ~ 0.2 GeV 2 



(A.3) 



(A.4) 



This matrix element was calculated independently and for different purposes in [K| 



and p7 |. On the other hand, one can introduce the two-particle DA's of twist 4 by 
the relation 



(0\d(0) W u(x)\7r(q)) = iU Qfl I due-^^+x'g^u)) 

-1 



~^~fiv I X [A X 



2 In 



qx / JO 



-luqx 



92(u) 



(A.5) 



where 4>(u) is the leading twist DA ( with (p as (u) = 6uu ) and all logs of x 2 are 
included in DA's. The twist 4 DA's gi(u),g2(u) can be now expressed in terms of 
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the three-particle DA's $y, The easiest way to do this is to use the equations of 
motion pi[ 

d(-x)^^u(x) = i J dvvd(-x)x a gG a/3 (vx)jpj 5 u(x) 



dx^ 



(d(-x)7 M 7 5 u(x)) = z y dvd(-x)x a gG a p(vx)^ 5 u(x) (A.6) 

(here <9 M stands for the translation operator). Taking the matrix elements and mak- 
ing use of the definitions ( A.lj) , ( |A.5| ) we obtain to this accuracy 

/ \ ^ r2 2 — 2 

9n u ) = 2 uu 

92{u) = y£ 2 ww(w- u) (A.7) 
G 2 {u) = h 2 u 2 u 2 

For further convenience we have introduced here the DA G2(u) defined by g2{u) = 
(d/du)G2(u). Now we are in position to show that for b 2 = —x 2 the combination 
b 2 uu is exactly the expansion parameter in the coordinate space. To see this, we 
choose the frame with x = (x + ,x_,x±) = (0,x_,x^). Then 

(0\d(0)-f +l5 u(0,x_,x ± )\n(q)) = if n q + f due~ m ^- (<j>(u) + x 2 ( 9l + G 2 )(u) 

Jo 



iUq+ f 1 due- iuq+x -6uu(l - ^-5 2 b 2 uu + ■ • •) (A.8) 
Jo 9 



It is seen that we have obtained the two first terms in the expansion of the pion WF 
^(^ = 2u — l,b 2 ) in powers of b 2 . It is now easy to find the value of {k 2 }^ from 
Eq.flOD : 

(0|J 7+75 (i V±) 2 u\rc(q)) = iUq + h 2 = if n q + (k 2 ± ) n , (^) 7r ~ (330 MeV) 2 (A.9) 

y 

We have thus demonstrated two important facts: 

1. The light-cone WF ^(u, k 2 ^) and higher twist DA's containing gluons explicitly 
are two sides of the same coin, as they are related by exact equations of motion. 

2. The combination k\/uu is a single variable which the WF depends on : 
^(u, k\) = ^(k 2 ± /uu). 

The same line of reasoning can be extended to the case of the p-meson. We therefore 
conjecture that the ansatz ([TJD is valid in the whole range of £ to the leading con- 
formal spin accuracy, i.e. it corresponds to a series of asymptotic two-particle DA's 
of higher twists. This approximation seems reasonable provided the WF \l/(£, k 2 ^) is 
evaluated at a high normalization point fi 2 ~ a few GeV 2 . On the other hand, it is 
such fi 2 that enters Eq. flU]) following from the duality arguments. 
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Appendix B 



The aim of this Appendix is to present the method |2(J which allows to estimate 
some high- dimensional condensates with a reasonable accuracy. The importance of 
this question (apart from the pure theoretical interest concerning a structure of the 
QCD vacuum) for the phenomenological purposes was emphasized in Section 3. It 
was observed that the large magnitude for the high dimensional condensates implies 
a strong momentum fluctuations in the transverse direction. Correspondingly, it has 
a remarkable impact on the pre-asymptotic behavior of the electroproduction. 



Let us start with the formulation of the idea exploited in this method. Consider 
some correlation function i J dxe iqx (T{Ji(x), J2(0)}) at large — q 2 . If the currents 
Ji, J 2 are chosen in such a way that in the chiral limit the perturbative contribution 
is zero, we end up with the leading contribution in the form ^ plus some nontrivial 
function, but with small numerical coefficient due to the loop suppression. Such a 
behavior, from the point of view of the dispersion relations (where each resonance |r) 
contributes to the correlation function as _ g2+m2 • (Q| Ji|r)(r| J2IO)) implies that in the 
chiral limit we can keep only the n meson contribution and, therefore, we make an 
identification (O) ~ (0\Ji\n) (n\J 2 \0) ■ An accuracy of this equation can be estimated 
by calculating corrections to the leading term which is ~ These corrections, 
which usually do not exceed 10 4- 20%, are responsible for the contribution of the 
higher states and proportional to ~ (0| Ji\r)(r\ J2IO). Those contributions have the 
same order of magnitude as the corrections in OPE, i.e. about 10 4- 20%. If it were 
not true, the matching of the right hand side and left hand side of the corresponding 
sum rules would be not possible. Therefore, if we knew the 7r-meson matrix elements 
(0| J'i 1 7r ) ( 7T I J2IO) independently (from PCAC, for example), we would estimate the 
condensate (O) with a high enough accuracy. 

We demonstrate the method by considering the following correlation function 
ijdxe iqx (T{J^(x), J 2 (0)})=<fT(g 2 ), J? = d^ l5 gaGu, J 2 = u l5 gaGd (B.l) 

where oG = cr^G®^ and a^ u = i/2{^/^ u — 7^7^) . The leading contribution is 
determined by the condensate of interest: 

T{q 2 ) = - 2 (q(gaG) 2 q). (B.2) 



At the same time, both ir meson matrix elements which enter this correlator 
are known- one of them (0 1^7^75 g aGu\n) can be expressed in terms of the matrix 
element (|A.4j) which we already know, and another one {7i\u'j5gaGd\0) can be found 
from PCAC. Finally we arrive at the following equation: 

(qg 2 {aGfq) = A5 2 ■ (q(gaG)q) ~ ^ g ° G ^ ~ m^qq), (B.3) 
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where we used the following it meson matrix element 



(0\dgG^ v u\*{q)) = iUq,5 2 ~ z/^i^M, S 2 ~ 0.2 GW 2 . (B.4) 



This matrix element was calculated in Ref.]2"7| (in terms of 5 2 )and in Ref.[fLU[| (in 
terms of mixed vacuum condensate). Numerically these results are in the perfect 
agreement to each other. 



As we mentioned, an accuracy of Eq. ([B.3Q can be estimated by calculating cor- 
rections to the leading term ( P-2| ). The most important loop correction can be 



calculated explicitly |20[ and it turns out to be 30 times smaller than the main 
contribution (P-2j ); therefore, the precision of our estimate ( P.3| ) for the condensate 
(qg 2 {crG) 2 q) is determined mainly by the accuracy of the 7r meson matrix elements 
which, we believe, is reasonably high. If we describe a deviation from the factoriza- 
tion prescription for the condensate (qg 2 (crG) 2 q) by introducing the parameter K 
as a measure of the non-factorizability (K = 1 if the factorization would work) , we 
get the following estimate from Eq. ( |B.3| ) for the parameter K: 

(qg 2 (aG) 2 q) = j(g 2 Gl v )(qq), if ~ 3 + 4. (B.5) 

which is the main result of these calculations. To be on the safe side we use K ~ 3 
in all our numerical estimates in the text. 
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